Fragmentation and the Bose-glass phase transition of the disordered 1-D Bose gas 
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We investigate the superfluid-insulator quantum phase transition in a disordered ID Bose gas in 
the mean field limit, by studying the probability distribution of the density. The superfluid phase is 
characterized by a vanishing probability to have zero density, whereas a nonzero probability marks 
the insulator phase. This relation is derived analytically, and confirmed by a numerical study. 
This fragmentation criterion is particularly suited for detecting the phase transition in experiments. 
When a harmonic trap is included, the transition to the insulating phase can be extracted from the 
statistics of the local density distribution. 
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A disordered low-dimensional non-interacting system 
is expected to always be in an insulating state [1]. Many- 
body interactions, in the case of Bose particles, may in- 
duce a phase transition to a superfluid state. This tran- 
sition is conventionally defined through the superfluid 
fraction or by the long-range behaviour of the one-body 
density matrix In ID at zero temperature, this lat- 
ter is expected to decay algebraically in the superfluid 
(SF) phase, and exponentially in the Bose-glass (BG) 
phase. The versatility and tunability of ultracold atomic 
systems have motivated, in recent years, the study of 
this fundamental phenomenon in a low-dimensional dis- 
ordered Bose gas for which however neither the 
superfluid fraction nor the one-body density matrix are 
easily accessible in experiments. Several different ways 
of characterizin g this p hase transition have thus recently 
been discussed llOMl4l . In spite of the many experimen- 
tal attempts |Bl Il5l4rn | however, a clear evidence of the 
superfluid-insulator transition in ID Bose gases has not 
yet been obtained. 

The notion of fragmentation has been widely used in 
relation to Bose-Einstein condensation [3, [13 ■ I n the 
case of a disordered Bose gas, fragmentation has been 
frequently evoked as a criterion for the transition from 
superfluid to Bose glass phase 0, fl9| . To our knowl- 
edge, however, a rigorous definition of fragmentation of 
the density profile, and a proof of its relation to the quan- 
tum phase of the gas, are still lacking. 

A quantity that naturally characterizes fragmentation 
is the probability distribution of the density (PDD), i.e. 
P(po), in the limit po —> 0. We define a state to be 
fragmented if P(pa 0) is nonzero. In this Letter we 
show that in the weakly interacting regime the density 
distribution fragments at the transition from superfluid 
to Bose glass. This result bears a clear experimental ad- 
vantage, as the density profile of a Bose gas is a quantity 
that can be easily investigated, within the spatial reso- 
lution of the experimental apparatus. To make a link to 
current experiments, a realistic situation of a gas con- 
fined in a harmonic trap is examined [20| . In this case, 
the transition between a quasi-condensate and an insula- 



tor - and a possible spatial separation between the two 
phases - can be unveiled through a local investigation of 
the statistical distribution of the density profile. 

To establish a link between the PDD and the phase 
transition, we study the superfluid fraction of the Bose 
gas. This quantity can be characterized by evaluating the 
response of the system to a velocity field, that is equiv- 
alent to imposing twisted boundary conditions |2l| . In 
fact, the superfluid fraction is proportional to the differ- 
ence of the energies in the moving frame, Eq, and in the 
rest frame, E , as [22] 
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Here O is the total phase twist, m is the mass of the 
bosons, L is the length of the system and N is the num- 
ber of particles. A mean field model of the weakly in- 
teracting ID Bose gas requires a description in terms 
of a a density, p = pa + 5p, and a phase operator, 9, 
as $f(r) ~ e l9 ( x ) p(x). The validity of the Bogoliubov 
prescription, (5p)/po -C 1, ensures that po contains the 
relevant information about the density distribution. The 
ground state density is the solution of a Gross-Pitaevskii 
equation (GPE) Q 



[-h 2 dl/{2m) + V{x) + gpo(x)] V 'po(x) = py/ Po (x) : 

(2) 

where p is the chemical potential, V(x) is the external 
potential and g is the interaction constant. From the 
Gross-Pitaevskii energy functional, computed at leading 
order in the phase twist, the energy difference is given by 
the kinetic term 
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where 9 = (9). Minimizing this energy with the con- 
straint J V# = shows that the total superfluid fraction 
is related to the harmonic average of the density [12J, [24 1 
as 
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The convergence of the integral in Eq. ([JJ is determined 
by the behavior of -P(po) m the limit po — > 0. If we 
express P(po — > 0) — Pq, then the condition to be in the 
superfluid phase is (3 > 0. A nonzero value of P(0) on 
the contrary, implies the insulator phase. 

We compute the PDD numerically, by solving Eq. @ 
on finite size systems with periodic boundary conditions. 
Configuration average has been adopted to increase the 
precision of the statistical sampling. This analysis is done 
for Gauss-distributed and Gauss-correlated disorder, de- 
scribed by 

(V(x)V(x')}=A 2 g e~ L ^t, (5) 

where A g is the disorder amplitude and r\ g is the spatial 
correlation length that introduces an additional energy 
scale, E c = . The third energy entering the prob- 
lem is the interaction energy, U = gNo/L, with No the 
number of bosons in the ground state. 

Previously we have characterized the phase boundary 
in independent ways, through the study of the superfluid 
fraction and the one-body density matrix [HI [l2| ■ Two 
limiting cases have been identified [HI]: a Thomas-Fermi 
regime (TF), where E c -C U, and a white-noise limit 
(WN), marked by E c 3> U. In these regimes, the phase 
boundary obeys power-law relations, A/E c = C(U / 'Ec) 1 , 
with 7 respectively equal to 1 and 3/4 (cfr Fig. H}. 

Fig. Ula) shows the PDD for fixed A g = 12.8 E c and 
increasing interaction U = 25.6 E c , 35.84 E c , 46.08 E c 
(all in the TF regime). From our previous analy- 
sis [Hi [13] these three values correspond to the BG, phase 
boundary and SF phases respectively. In the homoge- 
neous case P(po) is expected to have a single peak at 
the value Po/ph — F where pu is the constant solu- 
tion of the homogeneous problem. The inclusion of a 
small disorder [25J broadens this peak, but P(po) pre- 
serves a vanishing tail for po — > (solid blue curve in 
Fig. QJa)) in the superfluid phase. For increasing disor- 
der the weight of the low-density part becomes more im- 
portant (green dashed curve in Fig. HJa)) until the phase 
boundary is eventually crossed and the PDD develops a 
finite component in the limit po — > (red dot-dashed 
line in Fig. (Ha)). Fig. QJ6) shows a similar analysis 
carried out for the WN regime, for A g = 0.016 E c and 
U = 0.0032 E c , 0.0048 E c , 0.0064 E c . As for the TF 
case, from our previous study these three values lie in the 
BG, phase boundary and SF phases respectively. 

Comparing panels (a) and (b) in Fig. [1] it is clear 
that the PDD has different shapes in the TF and WN 
regimes. However, in both cases, the fragmentation al- 
lows to differentiate between SF and BG phases. The 
numerical analysis summarized in Fig. [1] confirms the 
criterion stemming from eq. (j4]), namely that the SF 
fraction is nonzero if and only if P(po — >• 0) = 0. We have 
checked that the same conclusions hold also for a speckle 
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FIG. 1: (color online) Probability distribution of the density 
in two different regimes, (a) Towards the TF regime, A 9 = 
12.8 E c for U = 25.6 E c (red dot-dashed), 35.84 E c (green 
dashed), 46 E c (blue solid). (6) Towards the WN limit, A g = 
0.016 E c for U = 0.0032 E c (red dot-dashed), 0.0048 E c 
(green dashed), 0.0064 E c (blue solid), pn is a constant value, 
solution of the homogeneous case. 

potential. Again, the choice of the potential strongly af- 
fects the shape of the PDD, but the limiting behaviour 
P(po 0) is only determined by the phase of the gas. 

P(po) can be evaluated analytically in certain regimes. 
In the TF regime, when the kinetic term is negligible, 
the density follows the external potential according to 
the TF approximation as 

p {x) = [p - V{x)]/g if V(x) < n, 
po(x) =0 if V{x) > p. (6) 

In this regime the distribution of po(x) reproduces the 
distribution of the potential at any finite value, with an 
additional finite contribution in zero, given by the sum of 
the regions where V(x) > p. A very similar feature can 
be indeed noticed in the insulating case of Fig. QJa) (dot- 
dashed line) where the PPD has a gaussian-like shape 
with a peaked contribution in 0. Following the fragmen- 
tation argument, this case is always insulating as it can 
be expected since the absence of a kinetic component 
prevents the formation of any quasi-long range order or 
superfluid flow. This is not in contradiction with the 
phase transition found for large values of U/E c , because 
the kinetic energy corrections to ((6|) are responsible for 
the build-up of the quasi-long-range order. 

Homogeneous systems are useful theoretical tools to 
inspect the phase transition, but current experiments 
[HI. Ilij are performed on trapped systems. The harmonic 
trap introduces a spatial inhomogeneity in the interaction 
energy U(x) — gpo(x). To deal with a realistic situation 
we consider a speckle disorder on the top of an harmonic 
trapping potential, as it has been investigated in recent 
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FIG. 2: Local probability distribution of the density in a har- 
monic trap (hut = 0.064 E c ) for three different values of 
the interaction energy Uc(0): (a) 3.2 E c , (b) 3.84 E c , (c) 
6.4 E c and for fixed speckle disorder intensity (A s = 3.2 E c ). 
The distributions are obtained from configuration averages. 
The red solid lines are the interaction energy profiles in the 
disorderless cases. Blue dashed lines represent the spatial dis- 
tribution of the average interaction energies. The green thin 
lines indicate the boundary in the thermodynamic limit. 
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and we perform the simulation on 
a system of experimentally achievable size [13] (200 cor- 
relation lengths). The speckle potential results from the 
interference pattern produced by the scattering of coher- 
ent light through a rough plate. Its amplitude, A s , is 
defined as the standard deviation of the potential and 
its distribution is bound from one side (above or below 
depending on the detuning of the laser) and decays expo- 
nentially on the opposite side [27J • Its correlation length 
•q s is given by the numerical aperture of the focusing lens, 
and in particular is linked to the cutoff in fc-space k c as 
% = l/k c . 

The contour plots in Fig. [2] show the spatially resolved 
probability distributions. P x (po), as a function of po- 
sition and interaction energy (U(x)) for fixed disorder 
amplitude, A s = 3.2 E c . The blue dashed line is the av- 
erage interaction energy U(x) = J gpo{x)P x (po)dpo, and 
the red solid line represents the interaction profile in ab- 
sence of disorder Uc(x), i.e. the solution of the GPE 
with V{x) = muj^x 2 /2, where uit is the trapping fre- 
quency. As we are considering a speckle potential bound 
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FIG. 3: Effect of a finite resolution on the local probabil- 
ity distribution of the density in a harmonic trap (hwt = 
0.064 E c ) for the insulating state of Fig. Ufa): U c {0) = 3.2 E c 
and A 3 = 3.2 E c . Panel (a) shows a resolution R — 1, whereas 
in panel (b) R = 2. 



from above (attractive), the distribution P x {po) is not 
centered around its mean value. Note also that the dis- 
order combined with a smooth trapping potential makes 
the average density profile different from the clean one, in 
fact it extends beyond the disorderless TF radius and it 
reaches a lower value in the center of the trap. This fact 
is relevant when comparing the critical values of interac- 
tion at the transition between the trapped and untrapped 
cases. 

It is important to understand whether the fragmenta- 
tion threshold, in the trapped case, coincides with the 
one extracted in the thermodynamic limit. For the cho- 
sen value of the disorder amplitude, the fragmentation 
analysis in the untrapped case gives a transition at an 
interaction strength Uh — 3.1 E c . In panel (a) the gas 
is fragmented at all positions, consistent with the ho- 
mogeneous boundary, as the maximal interaction energy 
£7(0) ~ 2.7 E c < Uh- When the interaction energy is 
increased, SF and BG phases coexist (see panels (6) and 
(c)). The fragmentation threshold of the interaction en- 
ergy can be extracted directly from the spatially resolved 
PDDs, these values for the average interaction at the 
threshold coincide with Uh (cfr. Fig. [^6) and (c)). This 
is true as long as the local density approximation is valid 
(ftwt <C U). For tighter traps (simulations not shown), 
the fragmentation line moves to lower values of the aver- 
age interaction energy, that can be interpreted as a pen- 
etration of the SF into the insulating state (cfr. Fig. @|. 
This is obvious when considering the limit tkOt 3> U, A, 
where the density profile is the non-fragmented harmonic 
oscillator ground state. The penetration effect is more 
pronounced when the transition occurs at the edges of 
the trap where the density gradient is steeper. There- 
fore, the best agreement between the phase transition of 
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the trapped and homogeneous gases is obtained when the 
transition occurs close to the center of the trap. 

In current experiments the resolution of the apparatus 
represents the main limitation to the correct interpreta- 
tion of the experimental data. For this purpose we study 
the role of a finite resolution on the statistical analysis 
that we propose in this article. The resolution is included 
in the simulation via a convolution of the ground state 
wavefunction with a gaussian of standard deviation R. 
The finite resolution makes it harder to identify the in- 
sulator phase through a statistical study of the density 
because it smooths out the profile, cutting narrow min- 
ima from the statistics. This limitation acts differently 
in different regimes: supposing the ratio r\jR fixed, the 
finite resolution would not affect the statistics in the WN 
limit, where modulations of the density occur on the scale 
of the healing length, much longer than the resolution. 
On the contrary, if r\jR < 1, this would strongly affect 
the statistics in the TF regime, where the lengthscale of 
the modulation would be the typical correlation length of 
the potential. Fig. [3] shows the effect of finite resolution 
on the probability distribution of a trapped Bose gas for 
U c {0) = 3.2 E c and A s = 3.2 E c . The statistical analy- 
sis of the distribution (Fig. [2ja)) identifies a completely 
fragmented state. The introduction of a finite resolution 
(R = 1) in Fig. G2a) opens a window in the center of 
the trap as if the Bose gas was partially superfluid. This 
effect become even more significant for a worse resolution 
(R = 2), as shown in Fig. [3jb), where the gas appears 
to be predominantly superfluid. For this set of param- 
eters, the transition would be placed at U ~ 3, 2.9 E c 
respectively. 

Both the finite resolution and the harmonic trap tend 
to displace the apparent fragmentation threshold towards 
the insulator phase with respect to the actual boundary. 
Therefore, in experiments, the superfluid fraction is over- 
estimated, but a quantitative correction to the boundary 
can be computed knowing the value of the resolution and 
the intensity of the trap. A qualitative sketch of the ef- 
fect of real conditions on the actual phase diagram is 
shown in Fig. |4] From this sketch it appears that the 
trapped system does no longer reflect the properties of 
the thermodynamic limit deep in the WN regime. In- 
deed, fragmentation would occur on a scale much longer 
than the ground state of the trap. 

In conclusion, we have demonstrated a relation be- 
tween the quantum phase transition of the ID Bose gas 
and the probability distribution of the density in the 
mean field limit. We found that the superfluid phase 
is marked by a vanishing probability at zero density 
whereas the insulating phase develops a finite component 
in the limit P(po — > 0). This gives a viable route to the 
determination of the phase of the gas, and therefore of 
the phase transition, through a statistical analysis of the 
density profiles. 

We thank B. Deissler, R. Hulet, P. Lugan, E. Runge 
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FIG. 4: (colors online) Phase diagram of the ID Bose gas 
in presence of disorder. The effect of realistic limitations in 
experiments are qualitatively shown. The black thick line 
represents the actual phase boundary, the blue dotted line 
is the phase boundary extracted with a finite resolution of 
the experimental setup, for a fixed ratio rj/R (see text). The 
red dashed line is the fragmentation boundary considering 
a harmonic potential of fixed trapping frequency. The trap 
tends to displace the boudary toward the insulating phase. 
For hu)t > U the local density approximation breaks down 
and the profile is never fragmented. The combined effect is 
shown by the green thin line. 
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